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ABSTRACT 

We study the abundance of satellites akin to the brightest, classical dwarf spheroidals around galaxies 
similar in magnitude and isolation to the Milky Way and M31 in the Sloan Digital Sky Survey. From a 
combination of photometric and spectroscopic redshifts, we bound the mean and the intrinsic scatter 
in the number of satellites down to ten magnitudes fainter than the Milky Way. Restricting to 
magnitudes brighter than Sagittarius, we show that the Milky Way is not a significant statistical 
outlier in its population of classical dwarf spheroidals. At fainter magnitudes, we find an upper limit 
of 13 on the mean number of satellites brighter than the Fornax dwarf spheroidal. Methods to improve 
these limits that utilize full photometric redshift distributions hold promise, but are currently limited 
by incompleteness at the very lowest redshifts. Theoretical models are left to explain why the majority 
of dark matter subhalos that orbit Milky Way-like galaxies are inefficient at making galaxies at the 
luminosity scale of the brightest dwarf spheroidals, or why these subhalos predicted by ACDM do not 
exist. 

Subject headings: dark matter — galaxies: dwarf — galaxies: formation — halo — Local Group 



1. INTRODUCTION 

Around the Milky Way (MW) Galaxy orbit nearly two 
dozen known satellite galaxies that have a range of mag- 
nitudes down to twenty times fainter than the MW it- 
self. The two brightest satellites, the Large Magellanic 
Cloud (LMC) and the SmaU Magellanic Cloud (SMC), 
are approximately two and four magnitudes fainter than 
the MW, respectively. Numerical simulations indicate 
that there is a < 10% chance that a IO^^Mq, MW- 
mass dark matter halo hos ts two satellites as massive 
as the Magellanic Clouds (jBovlan-Kolchin et al.l 120101 : 
iBusha et al.ll20lTh . From an observational perspective, 
analysis of galaxies from the Sloan Digital Sky Survey 
Seventh Data Release (SDSS DR7) indicates that MW 
magnitude-galaxies have two galaxies as bright as the 
Magellanic clouds only ~ 5% of the time, and on aver- 
age h ave ~ 0.3 satellites within four magnitudes of the 
MW (iLhTet al.ll201lHGuo et al.l[20Tlt iLares et al.ll2?nTI : 
iToUerud et al.ll2011| ). Thus the MW is somewhat atyp- 
ical in having two satellites with this luminosity differ- 
ence, but the probability of having a given number at 
this scale is in excellent agreement between theory an d 
observation (see further discussion in lBusha et al.|[2"01lD . 

Beyond the Magellanic Clouds (MCs), the brightest 
satellites are the Sagittarius, Fornax, and Leo I dwarf 
spheroidals (dSphs), approximately six, eight, and nine 
magnitudes fainter than the MW, respectively. Though 
at present challenging, a measurement of the abundance 
of satellite galaxies as bright as these classical dSphs is 
important for several reasons. From the observational 
perspective, it improves our understanding of the MW 
and its place on a cosmic scale, providing information 
on the number of bright satellites around other galax- 
ies over a regime in wh ich we are believe d to be com- 
plete around the MW (jWalsh et al.jl2009D. in particu- 
lar away from the Galactic plane ( Klevna et al.l[T997l : 
lWillmanll2010( ). From the theoretical perspective, under- 
standing the abundance of classical dSphs is important 
because these objects reside in the least massive dark 



matter halos that contain visible light, corresponding to 
halo mass scales at which various processes including su- 
pernova feedback effects and suppression of gas accre- 
tion from reionization be c omes important ([Bullock et al.l 
[2000HBmson et al.l[200llSomervillel[200l . 

Though the above processes certainly affect the for- 
mation of dSphs, theoretical models that include them 
have still found it challenging to match both the distri- 
bution of luminosities of t he bright MW s atellites and 
their kinematic properties (jFont et al.l 120 lit ). Estimates 
of the luminosity function down to scales of the classi- 
cal dSphs predict that they reside i n dark matter halo s 
with velocity dispersions > 20 km/s (jCooper et al.l20 1^. 
These models typically predict tens of satellites brighter 
than Fornax, ~ l(f Lq. Both of these predictions are 
in tension with observations of the MW satellite popula- 
tion. On the one hand, the predicted velocity dispersion 
for bright satellites is nearly two times larger than the 
observed ^ 10 km/s velocity dispersion of bright dSphs. 
On the other hand, since the observational sample is com- 
plete for objects with reasonable surface brightness with 
luminosity > lO'^L©, it is difficult to invoke that a pop- 
ulation of objects this luminous has been missed by ob- 
servations. 

More detailed analysis of dSph kinematics indicates 
that their maximum circular velocities are < 30 km/s 
(jStrigari et al.l [2010l ). In particular, even though it is 
very bright, Fornax has a strongly constrained maxi- 
mum circular velocity at a relatively low ~ 20 km/s. 
Matching these results with numerical simulations of 
the Galactic satellite population in a ACDM cosmol- 
ogy (jSpringel et al.ll2008t iDiemand et"alll2008| ) indicates 
that on average there should be ~ 25 — 75 of dark 
matter satellites with maximum circular velocity greater 
than that of Fornax, that are either too faint to be 
detect ed in surveys or devoid of baryonic material en- 
tirely (jBovlan-Kolchin et al.ir20lH) . 

The above results indicate that from a combination 
of observations and theory, the classical problem of the 
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abundance of satellites within ACDM can be boiled down 
to a "Fornax problem": more specifically, why is it 
that for an observed galaxy at the luminosity scale of 
~ IO^Lq, there are scores of dark subhalos that have 
the same dark matter mass, but apparently no stars in 
them at all? There are a few outstanding ideas that re- 
main to answer this question. First, the tension may 
point directly to severe inefficiency and stochasticity of 
galaxy formation at the dark matter halo mass scale 
of the dSphs. Second, it may be that there are re- 
quired modifications to the current sample of numer- 
ical simulations: in particular, baryons may signifi- 
cantly inodifY_tti£_da£^jnMtCT 

Ce.g. IWadepuhl fc Springell (IMl : IParrv et alJ (IMl : 
Idi Cintio et all (|201lb V or there may even be a neces- 
sary m odification to the cosmological model (jLovell et alJ 
120 lit ). The third idea is perhaps the most straightfor- 
ward of all; namely, that the MW is rare amongst galax- 
ies of its kind, a result of a rare downward fiuctuation in 
its population of bright satellites. 

In this paper, we address for the first time this lat- 
ter issue of the distribution of satellites as bright as the 
classical dSphs around MW-analog galaxies. We con- 
struct an observational sample of spectroscopic galaxies 
like the MW using SDSS, and search for faint satellites 
using SDSS DR8 imaging data and photometric redshift 
probability distributions. At the faintest end, we place 
an upper limit on the number of satellites down to the 
magnitude scale of the Leo I dSph. At the brightest end, 
we find an average of < 2 objects brighter than Sagittar- 
ius. Our results imply that, in terms of its bright dSph 
satellite population, the MW does not stand out as a 
significant statistical outlier. 

2. DATA 

Our analysis begins by selecting MW-analo g galaxies, 
which is similar to the analysis presented in Liu et aH 
(|2011| ). We refer to these galaxies as pri maries. We use 
a ma gnitude for the MW of My = —20.9 (jvan c en Bergh 
2000|) . which corresponds to Mr = -20.4 (see 
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for discussion), and we consider primaries within 
±0.25 magnitudes of this value. We consider a second set 
of slightly brighter galaxies, with Mr = —20.7, for com- 
parison with M31. Primary galaxies are selected from 
the NASA-Sloan Atlas (vO_l_lJ3, which combines SDSS 
data for nearby galaxies {z < 0.055) with more complete 
spectroscopy and imaging from GALEX. We limit our 
analysis to isolated galaxies that have no galaxy with 
a luminosity equal to or greater than the MW or M31 
wi thin 0.4 Mpc ( following the same methods described 
in lLiu et al.ll201ll ). We make a few cuts to exclude very 
nearby galaxies which may have incorrect photometry, 
specifically those with r < 11, SIZE < 5*SERSIC_TH50, 
or with g — r < 0. 

To perform the search for satellites around MW and 
M31-analogs, we utilize the both DR 8 imaging and spec- 
troscopic catalog (|Aihara et al.ll201ll ) . We use the galaxy 
catalog and t he photometric redshi fts probability distri- 
butions from iSheldon et al.l ()2011[ ). which is complete 
for r < 21.8. The spectroscopic sample is complete for 
r < 17.77. We exclude primary galaxies for which the 
DR8 photometry is incomplete within 250 kpc of the 

^ http://www.nsatlas.org/data 



primary (using a mask kindly provided by Erin Sheldon; 
the motivation for this cut will become more evident in 
Sec. [3]). 

For a given magnitude difference from the primary. 
Am, and satellite absolute magnitude threshold. Table [T] 
provides the corresponding apparent r-band magnitude 
threshold and the number of primary galaxies. Columns 
4-6 correspond to our photometric sample, while columns 
7-9 correspond to our spectroscopic sample. The mean- 
ing of these samples will become more clear from our 
discussion in Sec.|31 For the photometric sample of MW- 
analog primaries, we have explicitly included Am = 7.7, 
correspon ding to the Fornax dSph. In comparison to the 
sample of iLiu et al.l (1201 ID . our sample with a thresh- 
old cut at the magnitude of the brightest satellites, i.e. 
the MCs, is smaller because we restrict to more nearby 
galaxies in the NASA-Sloan Atlas (z < 0.055). Since the 
primary focus of our analysis is on the faintest satellites, 
our conculsions are insensitive to the smaller sample of 
very bright satellites. 

3. METHODS 

Consider a galaxy that has been selected as a MW 
or M31-analog via aforementioned methods. We define 
the "signal" region as circular area centered around the 
primary, corresponding to a physical radius R at the 
redshift of the primary. For the main analysis in this 
paper we take a R = 250 kpc for the signal region, cor- 
respondin g to the appro ximate viral radius of the MW 
and M31 (jSpringel et al, 2008: Diemand et al.ll2008D . In 
the signal region we search for galaxies that are between 
two magnitudes fainter and a threshold of Am magni- 
tudes fainter than the primary. For the i*'' primary, we 
label the number of galaxies in the signal region as nj.,. 
For comparison to the number of counts within the sig- 
nal region, wc associate with each primary a background 
region. The background region is an annulus with an 
inner radius R, and an outer radius chosen to enclose 
the same area as the signal region. We label the number 
of galaxies within the background region associated with 
the z*'' primary as rift j. A local estimation of the back- 
ground via an annulus connected to the signal region has 
been show n to be unbiased, even accounting for galaxy 
clustering (fChen et al.l l200a ILiu et al.|[20TTt iLares et"al] 
l20TlHGuoet al.ll2011h . 

We are interested in obtaining an estimate for the mean 
number of satellites around MW and M31-like galax- 
ies, given the measurements of nt.i and n^., around a 
large sample of primaries. We obtain three separate es- 
timates for the combination of Ut^i and Uh^t, which differ 
both in the cuts that are made on the data sample and 
the method in which the sample is obtained. The first 
estimate uses spectroscopic redshift information for the 
satellites, while the second and third methods use pho- 
tometric redshifts. We now detail specifically how each 
of these estimates are obtained in turn. 

3.1. Method 1: Spectroscopic satellites 

Our first method for estimating rit^i and uses 
spectroscopic redshifts for both primaries and satel- 
lites. Since the DR8 sample is spectroscopically-complete 
down to r < 17.77, we are able to estimate the number of 
satellites brighter than Am = 6, 5, 4. As indicated in Ta- 
ble [U these magnitude differences correspond to primary 
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Table 1 

Properties of primary galaxies for each magnitude limit. The upper eight rows correspond to Milky Way-like primaries, the bottom seven rows 

correspond to M31-Iike primaries. 



Photometric Spectroscopic 



Sample Am Satellite Absolute Primary Apparent Number of Mean Number Primary Apparent Number of Mean Number 
Magnitude Magnitude Primaries of Satellites Magnitude Primaries of Satellites 



MW 


10 


-10.9 


11.8 


49 


< 400 










9 


-11.9 


12.8 


148 


< 100 










8 


-12.9 


13.8 


727 


< 30 










7.7 


-13.2 


14.1 


1117 


< 13 










7 


-13.9 


14.8 


2727 


3.10 ± 1.60 










6 


-14.9 


15.8 


7574 


1.50 ±0.47 


11.8 


49 


1.42 ± 0.56 




5 


-15.9 


16.8 


8388 


0.92 ± 0.26 


12.8 


148 


0.62 ± 0.18 




4 


-16.9 


17.8 


8388 


0.51 ± 0.15 


13.8 


727 


0.25 ± 0.05 


M31 


10 


-11.2 


11.8 


41 


< 300 










9 


-12.2 


12.8 


151 


< 60 










8 


-13.2 


13.8 


834 


< 13 










7 


-14.2 


14.8 


3370 


2.41 ± 1.55 










6 


-15.2 


15.8 


7427 


1.31 ± 0.49 


11.8 


41 


1.01 ± 0.54 




5 


-16.2 


16.8 


7432 


1.02 ± 0.47 


12.8 


151 


0.37 ±0.12 




4 


-17.2 


17.8 


7432 


0.50 ±0.21 


13.8 


834 


0.20 ±0.04 



apparent magnitudes r < 11.8,12.8,13.8, respectively. 
Specifically we find 49, f48, and 727 primaries that sat- 
isfy r < 11.8, 12.8, f3.8, respectively. For a given set of 
primaries, we select galaxies in the spectroscopic sample 
that satisfy two criteria. First, as described above we 
determine those galaxies in the signal region that satisfy 
R < 250 kpc, as well as galaxies within the correspond- 
ing background annulus. Second, we impose a cut so that 
both the galaxies in the signal region and the galaxies in 
the background region lie within a redshift Az = 0.001 
(300kms-^) of the primary. This redshift cut is appro- 
priate when accounting for both the expected physical 
size of the dark matter halos of the primaries, and for 
redshift space distortions. 

Using spectroscopic redshifts to identify satellites 
clearly reduces background contamination from galax- 
ies at vastly different redshifts than the primary. How- 
ever, the obvious downside is that, when demanding the 
satellites have a spectroscopic redshift, the sample of pri- 
maries is much smaller than the corresponding sample of 
primaries obtained by only demanding complete photom- 
etry for the satellites. Further, when restricting to spec- 
troscopic satellites we run out of very bright primaries, 
meaning we are unable to probe very faint magnitude 
differences. Am > 6, that are the main goal of this work. 

3.2. Method 2: Full photometric redshift distributions 

For our second estimate oirit.i and n^.a , we us e the pho- 
tometric galaxy sample of lSheldon et al.l ()201lD , which is 
complete down to r < 21.8. In this case, we are able 
to estimate the number of satellites down to magnitudes 
as faint as Am = 10, corresponding to primaries with 
r < 11.8. The total number of primaries for each Am is 
shown under the photometric h eading in Table [J 

For each galaxy with r < 21.8. lSheldon et all (|2011[) as- 
sign a probability, p{z), that it resides in one of Nphot = 
35 redshift bins. Method 2 specifically uses the informa- 
tion in these p{z) distributions. With these distributions 
in hand, again we start by considering a single primary 
galaxy. We use the same redshift bins in which the p(2:)'s 
are calculated, and from these locate the redshift bin that 
contains the primary. As described above we again search 
for galaxies within the signal region R < 250 kpc. Given 



this set of galaxies in the signal region, we then sample 
each of their corresponding discrete p{z) distributions. 
This provides us with a redshift for each galaxy within 
the signal region. From a single sampling of each p(z) 
distribution, we determine the total number of galaxies 
that fall into the redshift bin that contains the primary. 
To obtain a statistically robust estimate of the total sig- 
nal and background probability distributions, we repeat 
the sampling procedure ^ 100 times to obtain the mean 
number of counts in the signal region in the same redshift 
bin as the primary, {nt^i). Then following a similar rea- 
soning as above, we repeat this procedure for the annu- 
lar background region to obtain the corresponding mean 
number of counts in the redshift bin associated with the 
background for each primary, (nt^j). 

As long as the photometric redshift distribution for 
each galaxy is unbiased, the estimate of the number of 
satellites obtained from this method is expected to be 
just as accurate as the estimate that uses spectroscopic 
redshifts. However, systematics in photometric redshift 
information may bias the results. Of parti cular relevance 
is the fact that the training set used in I Sheldon et al.l 
(|2011t ) to obtain the photometric redshifts may not have 
a fully representative sample of faint, low redshift galax- 
ies. This is precisely the sample of galaxies that we are in- 
terested in identifying as satellites. If the training sample 
is underrepresented at low redshift, this will likely lead 
to a biased p{z) estimate for potential satellites, and un- 
derestimate of the satellite number associated with low 
redshift hosts. As shown below, a comparison of the re- 
sults obtained from this method to those obtained from 
the spectroscopic sample (method 1), over the regimes 
in which both are complete, may provide an estimate of 
the magnitude regime in which photometric redshifts are 
not representative for faint, low redshifts galaxies. 

3.3. Method 3: Mean photometric redshifts 

Our third method to estimate nt^i and n},^^ also uses 
information contained in the photometric redshift distri- 
butions. However, rather than directly sampling the full 
p[z) distributions as above, this method simply uses the 
information contained in the mean of these distributions. 
Once a mean photometric redshift is assigned to a galaxy, 
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Figure 1. Mean number of satellites within Am of MW-like primaries, using the three different methods described in Sec.|3] In all panels, 
blue diamonds include only satellites with spectroscopic redshift within Az = 0.001 of the primary (method 1), and red triangles sample 
the full photometric redshift probability distributions (method 2). In the three left panels, black squares only include satellites with mean 
photometric redshift less than 0.2, while in the right-most panel black squares only include satellites with mean photometric redshift less 
than 0.5 (method 3). In each of the four panels, the magnitude limit on the population of satellites from the photometric analyses (methods 
2 and 3) are indicated. In all panels, upper limits are indicated as downward arrows. 



implementation of this analysis method boils down to 
determining an appropriate redshift cut on galaxies in- 
cluded in the signal and background region. 

To motivate the appropriate redshift cut. define the 
mean photometric redshift for a galaxy as 



Nphot 

E 



(1) 



where p{z.i) is the probability for the z*'* redshift bin. 
Though it is simple, ascribing to a galaxy a single redshift 
based upon Eq.[l]will typically overestimate the true red- 
shift of the galaxy. More specifically we are interested in 
the mean photometric redshifts for satellites that span 
the redshift range of our primary galaxies; this corre- 
sponding redshift range is 0.055. Examining the faintest 
sample of training set galaxies, 20.8 < r < 21.8, over this 
redshift range of the primaries, we find that the ratio of 
the mean photometric redshift to the true redshift, ztrue, 
spans a wide range {zp) / ztme ~ 3 — 30. Further, even for 
the brightest sample of training set galaxies fainter than 
the spectroscopic completeness limits, 17.8 < r < 18.8, 
we find {zp) / ztme ~ 1 - 10. 

In spite of the fact that typically {zp) / ztme ^ 1, it is 
still possible to use the information on the mean photo- 
metric redshifts, provided that we place an appropriate 
upper cut on the mean photometric redshift of a galaxy 
that is allowed in our sample. In other words, we must 
estimate the number of true satellites that are lost from 
the sample when including galaxies less than a given {zp). 

We obtain this estimate for the loss fraction by again 
considering training set galaxies over the redshift range 
spanned by the set of primary galaxies of interest. Specif- 
ically, out to the maximum primary redshift, we examine 
the distribution of {zp) for training set galaxies that are 
brighter than a given apparent magnitude. Then from 
this distribution we determine the value of (z. 



■p/90: 



which 



we define as the redshift below which 90% of the galaxies 
reside. We determine (zp)9o for each apparent magnitude 
threshold cut. For example we find that when searching 



for satellites down to r < 21.8, {zp)'^^ — 0.5. Further, 
we find that for r < 19.8, {zp)go = 0.2. These are the 
cuts that we utilize in our analysis below. Because of the 
large sample of primaries available when using the pho- 
tometric catalogue, we find that the results of method 3 
are generally insensitive to the chosen value for (2^)90. 

3.4. Likelihood 

With the methods outlined for determining the data 
sample of galaxies in the signal and background regions, 
it remains to construct the likelihood for the data. For a 
total sample of Hp primaries, the mean number of satel- 
lites is the difference between the mean of the signal and 
background counts. 



^ Tip _^ Tip 

= — Yl — E 



(2) 



1=1 



'P ,=1 



We define the variance of the signal and background 
distributions as and a^. Assuming that the num- 
ber of galaxies around a primary and the background 



are uncorrelated, the variance is crj 



(7? . Defin- 



ing the parameter set that we estimate from the data as 
X — [/iis, as, fib, CTb], for a given magnitude threshold Am, 
the probability for the mean and variance is 



rip 



P(.'|Am)ocn^^exp 



1=1 

X 



exp 



2al 



(3) 



The probability distribution P{x\/^m) can be thought of 
as the probability for the number of counts in the signal 
region, the first term in Eq. [3l weighted by a prior given 
by the number of background counts, the second term in 
Eq.H 

Equation[3]is a general formula that can be used to esti- 
mate the mean and the variance of the satellite probabil- 
ity distribution, independent of the number of primaries 
that contribute to the sample. We integrate Eq. [3] and 
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marginalize over the background mean and variance to 
obtain the probability for jig and as- In our results below, 
error bars are derived as the area centered on the mean 
containing 68% of the cumulative probability distribu- 
tion for /is. We use uniform priors on the parameters 
that we estimate, P{x) — const. 

Assuming that the respective number counts nt^i and 
rib^j for a given primary reflect the mean of an underlying 
redshift probability distribution, Eq. [3] follows directly 
from the central limit theorem. This is explicitly true 
in the case of method 2, which makes the replacements 
nt^i — > {^1,1) and rib^i — {ub^i) in Eq. |31 In the limit of 
a large number of primaries, we have verified that Eq. |3] 
reproduces the true mean of t he full satellite p robability 
distribution, as determined in lLiu et al.l (|2011[ ). 

4. RESULTS 

We begin by comparing the results of the various meth- 
ods for estimating the number of satellites for several 
samples of host galaxies. The four panels in Figure [1] 
show the mean number of satellites around MW-like pri- 
maries that are brighter than Am magnitudes fainter 
than the primary, /is(< Am), for the three methods de- 
scribed above. In each of the four panels, there is a dif- 
ferent threshold cut on the magnitude of signal and back- 
ground galaxies for the two methods that use photomet- 
ric redshifts (methods 2 and 3). From left to right, the 
magnitude cuts on signal and background galaxies are 
r < 17.8,18.8,19.8,21.8. Results from the spectroscopic 
analysis (method 1) are shown as blue diamonds in all 
panels; by definition this method only includes galaxies 
in the signal and background regions with r < 17.77. Re- 
sults from the full p{z) sampling analysis (method 2) are 
shown as red triangles, and from the cut on mean photo- 
metric redshifts (method 3) are shown as black squares. 

Figurc[T]indicatcs that for satellites with r < 17.8 (left- 
most panel), we find good agreement between all three 
methods, with large error bars for the method using sim- 
ple p{z) cuts. For these primaries, the p{z) distribution 
is determined by galaxies in the main sample of SDSS, so 
the training sample is expected to be fully representative. 
The agreement between methods 1 and 2 is maintained 
for Am = 4 and r < 18.8 (left two panels). However, we 
find that this agreement weakens for Am = 5,6, where 
the mean abundances as determined from method 2 un- 
derestimate the results from the spectroscopic analysis. 
This underestimation is most evident in the right panel of 
Fig. [U and as discussed above most likely indicates that 
galaxies with r > 18.8 and z < 0.055 are not adequately 
represented in the training sets. In addition, as expected, 
method 3 tends to overestimate the mean abundance for 
all Am and all magnitude limits. This is true both in the 
regime in which (zp)go = 0.2, for satellites with r < 19.8 
(left three panels), and {zp)go = 0.5, and for satellites 
with r < 21.8 (right panel). These results indicate that 
the photometric redshift distributions for satellites dim- 
mer than r = 17.8 are likely biased towards higher red- 
shifts, as a result of increasing incompleteness in the very 
lowest redshift galaxies in the training sample for dimmer 
galaxies. 

What leads to this bias in the photometric redshifts, 
and is it avoidable? First, it is important to realize that 
our regime of interest is very low redshift, z < 0.02, well 
outside the redshift range where photometric redshifts 



arc general ly used and tested. Th e p{z) distributions cal- 
culated bjJSheldo^e^n ()201lD use t raining sets from 
SDSS (lAihara et al.ll2011[). PR IMUS (ICoil et al.ll201lD 
zCOS MOS (ILiUv et al.l 120071) . 2SLAQ (ICannon et ah 



2np), VVDS (IGarilh et al.ll2008l), DEEP2 I Weiner et al 



CN0C2 (lYee et al.1 [20001) . C FRS (ILillv et al 



199^, and TKRS ijWirth et all |2003) . For galaxies in 



our main regime of interest, 17.8 < r < 21.8, the train- 
ing sample is dominated by PRIMUS (|Coil et al.ll201lL 
, Cool et al, in preparation), which covers more than 9 
sq. degrees to a depth of zab ^ 23.5 (5.2 sq. degrees 
arc included in the current training set). For PRIMUS, 
which is focused on science at z > 0.2, there arc at 
least two distinct issues that impact our analysis. First, 
the 4000 Abreak falls out of the wavelength regime for 
z < 0.2; without this feature the low- redshift spectra 
were more likely to be assigned a lower confidence flag 
{Q = 3) and were excluded in the p(z) estimation (Car- 
los Cunha, private communication). Second, the analysis 
pipeline does not even attempt to measure accurate red- 
shifts in the regime of our primaries (z < 0.02) because 
of the limited velocity resolution of the low-dispersion 
prism used by PRIMUS (John Moustakas & Alison Coil, 
private communication). Both of these will impact the 
completeness of the low redshift sample and are likely to 
bias the p{z) distribution against the redshift range of 
our primary samples. 

One might also consider whether the biased p{z) distri- 
bution could be avoided with another photometric red- 
shift code. However, this range is challenging for any 
algorithm. For example, template-fitting codes without 
priors have known failur e modes at very low redshift 
(jMandelbaum et al.ll2008| ) , and this would likely contam- 
inate the low redshift sample. 

From the information in Figure [1] we are able to obtain 
our best estimate for /is(< Am) over the entire range 
of Am. Specifically for our best estimate of fJ.s{< Am) 
we use the spectroscopic results for Am < 6, and the 
results from method 3 for Am > 6. Figure [2] shows these 
as the main results of our analysis, in comparison to the 
observational points from the MW and M31. The results 
from method 3 are explicitly shown in the 6th column of 
Table [Tl and the results from method 1 are shown in 
the last column of Table [TJ Note that for the cases with 
a small number of primaries, which for our MW sample 
corresponds to Am > 7, only an upper limit on the mean 
ca n be determ i ned. 

iLares et al] ()20lH ) report that the satellite population 
is uncertain in the central 100 kpc around primary galax- 
ies due to contamination from features in the extended 
halos of primaries. We have directly examined the images 
of many nearby, bright primaries and fi nd that gener- 
ally th is systematic is accounted for in the lSheldon et all 
(|2011f ) galaxy sample used in our analysis. In addition 
we have explicitly examined how our results are affected 
when excluding the 100 kpc around primaries, and we 
generally find minimal changes relative to those results 
presented in Fig. [2] For example, for Fornax-magnitude 
satellites around MW-like primaries, we find that the 
90% c.l. upper limit on the mean reduces from 13 to 9 
when excluding the central 100 kpc. 

In addition to the distribution of the mean number of 
satellites, from our likelihood analysis we are able to esti- 
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Figure 2. Left: Mean number of satellites brighter than Am magnitudes fainter than the primary galaxy, assuming primaries within 
±0.25 magnitudes of the Milky Way. Blue diamonds are determined from the spectroscopic sample of satellites (method 1), black squares 
from the photometric sample (method 3). The solid errors are the uncertainty on the mean, the thin, dashed errors are the intrinsic scatter 
{as from Eq. IS}. The arrows indicate 90% c.l. upper limits. The red triangles indicate the Milky Way satellites. Right: Same as left, except 
for primaries within ±0.25 magnitudes of M31. 



mate the intrinsic scatter, or as, for each magnitude bin. 
As for the mean, as is determined from the full probabil- 
ity density distribution using Eq. [3] As an example for 
MW-like primaries and Ajti = [4, 5], we find a mean in- 
trinsic scatter of cr^ = [0.56±0.04, 0.89±0.19], where the 
errors represent one-sigma uncertainties as above. The 
best-fitting values for as are shown as thin, dashed er- 
ror bars in Fig. for A m < 7. Via the method out- 
lined in iLiu et alT ()2011[ ) , we are also able to estimate 
the full probability distribution down to Am = 5; here 
we find that the probability to obtain [0,1,2,3,4] satel- 
lites with Am < 5 is [0.59,0.25,0.11,0.03,0.02]. Down 
to fainter magnitudes, the spectroscopic sample is too 
sparse to measure the full satellite probability distribu- 
tion. These results indicate that there is still substantial 
intrinsic scatter in the satellite population, even at the 
brightest scales. 

We note that the limits we present are strictly valid 
over the regime of surface brightness where the SDSS 
DR8 data is complete. For galaxies with half-light sur- 
face brightness < 22.5 mag/acrmin^, which is the sur- 
face brightness of Le o I, the SPSS spect roscopic data 
is ~ 90% complete (jBlanton et al.ll2005D . For galax- 
ies with surface brightness similar to Fornax or Sculptor, 
< 23.5 mag/acrmin^, the surface brightness complete- 
ness is ^ 50%. We note that the surface brightness 
incompleteness has not been estimated directly for the 
photometric sample, but the agreement our results be- 
tween the two sets of samples indicates that the issues 
may be of similar magnitude. We can get an estimate as 
to how the surface brightness incompleteness affects our 
resul ts by comp aring the integrated luminosity function 
of Blanton et al . (2005) that is corrected for incomplete- 
ness, as compared to the measured luminosity function. 
Down to the magnitude of Fornax, for example we find 
that the luminosity functions differ by a factor < 2 going 



to down surface brightness of 24 mag/acrmin^. If inter- 
ested in constraining the population of objects down to 
this surface brightness, this factor should be taken as a 
conservative correction to the limits that are presented 
in Fig. [21 Thus in order to obtain many more bright 
satellites than are observed in the MW, it is clear that 
these satellites must have surface brightness much dim- 
mer than the known bright MW satellites. 

5. COMPARISON TO PREVIOUS RESULTS 

There have been several recent analyses on the pop- 
ulation of bright satellites around MW-analog galaxies 
along the lines presented in this paper. It is instructive 
to compare the results presented here to these previous 
a nalyses. 

iGuo et al.l (|2011[ ) used SDSS DR7 to construct the lu- 
minosity function of satellites down to the magnitude 
scale of Fornax, correcting for the incompleteness of 
SDSS. These authors used best-fitting photometric red- 
shifts from DR7 to eliminate obvious background galax- 
ies. Our analysis differs from these authors in that we 
utilize both DR8 imaging and a maximum likelihood 
method that incorporates full photometric rcdshift prob- 
ability distributions. Wc also directly quantify the bias 
in abundance counts for faint satellites that is incurred 
when utilizing available photometric redshifts. Via some- 
what different me thods for cutting background galaxies, 
iLares et "all (|2011[ ) use DR7 data to obtain a mean num- 
ber of satellites down to the magnitude of Sagittarius for 
projected radii > 100 kpc. As we discuss above, we have 
verified that our results are consistent with these authors 
over the radial range considered, and further that we do 
not incur a significant bias by including g a laxies within 
projected radii < 100 kpc. iToUerud et al'l (|2011[ ) utilize 
the DR7 volume-limited spectroscopic sample and find 
that ~ 40% of MW-analogs have satellites brighter than 
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the LMC within 250 kpc. iJames fc Ivorvl ()2011[ ) use Ha 
narrow band imaging to search for start forming galax- 
ies around 143 spiral galaxies like the MW, and find that 
nearly two-thirds do not have satellites that resemble the 
Magellanic Clouds. These latter two results are consis- 
tent with the spectroscopic results that we present for 
bright satellites. 

6. DISCUSSION AND CONCLUSION 

We have used DR8 photometric redshift data to limit 
the mean number of satellites around MW-analog galax- 
ies down to ten magnitudes fainter than the MW. At 
least down to the scale of Sagittarius, the results indi- 
cate that the MW is not a significant statistical outlier 
in its number of bright, classical satellites. 

Our 90% c.l. upper bound of < 13 satellites brighter 
than the Fornax dSph already places a strict bound on 
the efficiency of galaxy formation at the dSph lumi- 
nosity scale. This is particularly true considering that 
there are anywhere from ~ 25 — 75 dark matter sub halos 
in the Aquarius simulations (jSpringel et al.l 120081 ) that 
have present-day circular velocities greater than that 
of Fornax. Surface brightness incompleteness could in- 
crease this number, but likely not enough to bring it 
into agreement with predictions for the number of dense 
satellites in simulations. However, it is very interest- 
ing to note that the observational result we present is 
perfectly consistent with abundance matching extrapola- 
tions for the satellite luminosity function, which predict 
~ 1.2, 1.7 satellites for magnitude differences Am = 7, 10 
(|Busha et al.ll20li[ ). This does not guarentee that such 
models will have the correct velocity function; in fact 
it appears increasingly difficult to simultaneously match 
both the luminosities and velocities of all of the satellites 
down to the Fornax scale. 

In the future it will be exceedingly important to in- 
crease the sample of primary galaxies around which it is 
possible to measure satellites as faint as Fornax. Measur- 
ing the magnitude distribution at this faint scale will go 
a long way towards determining if the mapping between 
bright dSphs and dark matter subhalos is revealing the 
presence of detailed baryoni c physics not yet accounted 
for in numerical simulations (iParrv et al.ll2011[) o r about 
the properties of dark matter ( Lovell et al.ll201ll ). 

The principle uncertainty in accurately determining 
the satellite distribution from the SDSS photometric 
sample is the fidelity of the photometric redshifts. We 
have shown that our method using the full photomet- 
ric redshift distribution is in excellent agreement with 
a method that directly uses spectroscopic redshifts for 
bright satellites where the training sample is representa- 
tive. However, this method is biased to lower satellite 
numbers for dimmer magnitudes, which is an indication 
that the spectroscopic training samples used to construct 
the photometric redshift distributions are systematically 
missing the lowest redshift galaxies. Our measurements 
would be substantially improved if these samples were 
unbiased, as this would allow us to use the p{z) method 
for a substantially larger number of host galaxies. It may 
be possible in the very near future to extend the train- 
ing sample to include data from GAMAQ, extending to 
r < 19.8, which would allow us to extend our measure- 



ment using the p{z) method to the Fornax scale. With 
some care, it may also be possible to refine and improve 
the PRIMUS redshift finder in the range relevant to our 
low-rcdshift primaries (z < 0.02), thereby providing a 
more reliable photometric redshift training sample. 

Given forthcoming data sets it will also be possible to 
significantly increase the number of primary galaxies; for 
example, the Dark Energy Survey (DESj3 will produce 
a survey over 5000 sq. degrees and will observe galaxies 
to 24.3 in the i-hand. If one considers satellites down 
to r = 24, the DES is expected to identify more than 
4000 primaries with satellites to Fornax magnitude dif- 
ferences, about 1600 primaries with Leo I-like satellites, 
and nearly 100 primaries with satellites as dim relative 
to their primary as Sculptor is to the Milky Way. At 
the Fornax scale, the statistics should be large enough to 
get a solid measurement using background subtraction 
even without robust photometric redshifts. The anal- 
ysis here indicates that the primary challenge for the 
dimmest satellites around these primaries will be deter- 
mining their redshift distribution. 

In the future it will also be important to obtain kine- 
matics on spectroscopically-confirmcd satellites, as well 
as those that have a high probability to be satellites from 
their photometric redshifts. Comparing the velocity dis- 
persion of these satellites to the ~ 10 km/s velocity dis- 
persions of the MW dSphs will allow for a determination 
of both the luminosity function and the mass function 
of satellites down to the scale of classical dSphs. This is 
only currently possible with the sample of MW dSphs. 
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